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Abstract 

A new string theory is proposed as a candidate for the large N Hmit 
of QCD. In this theory, strings are constrained to be non-backtracking; 
a condition that is essential for the gauge invariance of the underlying 
field theory. This condition is implemented by first placing the theory on 
a lattice and then introducing fermionic variables on the world sheet. The 
naive continuum limit leads to a generalized Thirring model on the world 
sheet, and it is suggested that the application of the renormalization group 
should drive the coupling constants to a conformally invariant fixed point. 
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1. Introduction 

It has been known for some time that in the large N hmit,QCD without mat- 
ter fields is equivalent to a free string theory(l). Despite considerable effort(2), 
this string theory is yet to be discovered. Recently,there has been renewed in- 
terest in this problem(3),as well as in other approaches to QCD at large N in 
general(4). Whatever it is, the QCD string cannot be the same as the standard 
Nambu-Polyakov string for at least the following reasons: 

a)There are no zero mass particles(and of course no tachyon)in the QCD spec- 
trum. 

b)The Regge trajectories of QCD should be curved, as opposed to the straight 
trajectories of the standard string models. This is true of even two dimensional 
QCD coupled to matter(5). 

One possible approach to this problem is to reformulate the dynamical 
equations of QCD in terms of gauge invariant quantities such as the Wilson 
loops. These equations have been derived for both the lattice(6) and the con- 
tinuum(7) versions of the theory. Unfortunately,this approach has not thus far 
led to a simple string picture. One of the problems with the lattice appraoch is 
that the resulting string picture is extremely complicated(8).This is probably 
due to the simplicity of the starting point, which is the usual Wilson single pla- 
quette interaction;it is conceivable that a more complicated interaction, which 
is in the same universality class and which therefore has the same continuum 
limit, may give rise to a simpler string. However, since it is not clear a priori how 
to determine the interaction that corresponds to a simple string picture on the 
lattice, it is tempting to bypass the lattice altogether and work directly in the 
continuum with Migdal-Makeenko equations(7).In this case, one encounters the 
problem of finding a suitable non-perturbative regulator to avoid divergences. So 
far, not much progress has been made along these lines. 

In this paper, we bypass the difficult problem of showing directly the equiv- 
alence between QCD and the corresponding string theory;instead,we focus on 
the constraints satisfied by Wilson loops. One can think of the world sheet of the 
QCD string as being traced by a spacelike Wilson loop propagating in time, and 
in the Hamiltonian description of this process, the physical Hilbert space is gen- 
erated by the collection of all spacelike Wilson loops. This forms an overcomplete 



set of states as a consequence of two kinds of constraints satisfied by the Wilson 
loops: 

a)The Mandelstam identities(9).They follow from the fact that N,the dimension 
of the color space ,is finite. It is generally assumed that in the limit N -^ oo,there 
are no identities of this type. Since we are always in the large N limit in this pa- 
per,we will not consider them any further. 

b) Every backtracking Wilson loop is equivalent to a non-backtracking one with 
the backtracking portion removed(see Fig. 1). We shall argue in the next section 
that this equivalence is crucial for the existence of an underlying gauge the- 
ory;if this condition is not satisfied, Wilson loops cannot be constructed from 
a gauge connection. Conversely,if this condition holds, there always exists a con- 
nection valued in some (in general infinite dimensional) color space. This is the 
condition we focus on here. 

In the standard bosonic string theory,one has an unrestricted sum over all 
possible world sheets,and therefore, taking a fixed time slice, all possible Wil- 
son loops, backtracking as well as non-backtracking, are present. Since the back- 
tracking loops are equivalent to the corresponding non-backtracking ones , there 
is clearly some overcounting.At the minimum, a string theory based on QCD 
should avoid this overcounting and sum only on world sheets whose cross sec- 
tions are non-backtracking loops. In this paper, we will study precisely such a 
string theory where the sum over the world sheets is suitably restricted. Because 
such a theory always admits a connection(see (b)above),it corresponds to some 
kind of gauge theory,and the question is, does it have the same interaction as 
QCD?In this paper,we shall leave this question open, and study this new string 
theory as an interesting theory in its own right. Since the underlying field theory 
is gauge invariant, it is a plausible candidate for the QCD string,hence the title 
of this paper. In a future publication(10),we hope to show that the version of the 
model where time is continuous and space is latticized corresponds to a gauge 
theory on a space lattice. However, the interaction is not the standard single pla- 
quette interaction of Wilson, but instead it is a suitably weighed infinite sum 
over the traces of all possible Wilson loops. This sum has the standard QCD 
interaction as its naive continuum limit, and if we believe this naive limit, which 
we afterall do in the case of the single plaquette interaction, the equivalence to 
QCD follows. 



The technical problem that will occupy us for most of the rest of this paper is 
how to impose the condition that strings are non-backtracking. Since it is difficult 
to handle this condition directly in the continuum,we first start with a regular 
space-time lattice and only at the end take the continuum limit. As a warmup ex- 
ercise, in section 2, we consider the simpler problem of non-backtracking random 
walk(Brownian motion), which involves summing over non-bactracking paths on 
the lattice. This problem can be solved in a number of different ways;the method 
which is easiest to generalize to sums over surfaces is the standard path integral 
approach. In order to implement the non-backtracking condition,we introduce 
additional fermionic variables, and we construct the action, making use of the 
transfer matrix method. The path integral must be supplemented by a boundary 
condition on the conserved fermion number to eliminate the unwanted back- 
tracking paths. In section 3, this method is generalized to non-backtracking sur- 
faces,which are defined as the direct product of two non-backtracking curves, and 
again a path integral formulation is obtained. The naive continuum limit of this 
model is discussed in section 4;the resulting two dimensional field theory that 
lives on the world sheet is a generalized Thirring model with the following prop- 
erties: 

a) It is scale invariant at the classical level. 
b)It is Lorentz invariant in the target space. 

c)It is Lorentz and parity invariant on the two dimensional world sheet as well. 
d)There is invariance under both vector and axial phase transformations of 
fermions, assuming absence of anomalies. As a result, there exists a chirally con- 
served fermion current on the world sheet. Just as in the case of random walk,suitable 
boundary conditions must be imposed on the conserved fermion current in order 
to eliminate the unwanted sectors of the model. 

There are several coupling constants in world sheet field theory that are 
easily calculated in the naive continuum limit. Of course, this is not the final 
answer;these constants should be allowed to run under the renormalization 
group equations and reach a fixed point. We are unable to carry out this anal- 
ysis,instead,we conjecture that there is a non-trivial fixed point and that the 
theory is conformally invariant at that point. Naturally,in searching for this 
fixed point, only those renormalization fiows that preserve the invariances a) 
to d) should be considered. The main result of this paper is therefore the fol- 



lowing: The original problem of non-backtracking strings has been transformed 
into the problem of finding the conformal fixed point of a generalized Thirring 
model invariant under the symmetries listed above. We feel that this is a result of 
interest, although finding such a fixed point, and solving the resulting model, are 
difficult problems, and they will not be attempted here. So far, in problems of 
this type, solutions have been found only in very special cases(ll), and hope- 
fully,methods developed in the study of conformal theories should prove useful 
in tackling this problem in the future. More advanced problems such as the value 
of the central charge and the higher order corrections in 1/A^ are also beyond 
the scope of this paper. 

2. Non-Backtracking Random Walk 

We start by fixing our notation and defining some terminology. Space-time 
(target space)is taken to be D dimensional;the value of D will not be fixed in 
this paper. The metric in the target space is fiat and can be either Euclidean 
or Minkowski. We are interested in both open and closed strings;we identify the 
open string with the Wilson line integral 

Ur{x,y) = P{exp r dx'.A{x')}, (2.1a) 

Jx 

where F is a path with endpoints x and y,A is the connection valued in some 
(infinite dimensional) matrix space, and P denotes path ordering.!! could be 
thought of as generating a line of fiux connecting two external quarks at x and 
y along the path F. Closed strings correspond to 

W = Tr{Ur{x,x)} (2.16) 

for a closed path F.The equivalence between the U's associated with the back- 
tracking path F and the corresponding non-backtracking path F pictured in 
Fig.l follows from 

Ur{x,y) = U-\{y,x), (2.2) 

where F~^ is F traversed in the opposite sense. As a result, the contributions to U 
from the segments a and b in Fig.l cancel. Of course, (2.2) follows immediately 
from (2.1a), but the converse is also true: If eq. (2.2) holds for arbitrary paths, 
and in addition, we also require the semi-group property 

f/r(i,2) = f/r(i)f^r(2) (2.3) 



for two arbitrary paths T^-^^ and T^'^^ that join at point y and give r*^-*^^) (see 
Fig. 2), then U can be written in terms in terms of a connection A as in (2.1). 
Therefore,eq's (2.2)and (2.3) guarantee the existence of an underlying gauge 
theory. In a random walk(or a string) problem,one has to deal with the paths 
directly, instead of the U's attached to them,and so it is convenient to define an 
abstract group generated by path multiplication (12). The group axioms are 

T{x,y)T-\y,x) = l (2.4a) 

and 

T^'^'\x,z)=T^'\x,y)T^'\y,z), (2.46) 

where, again, multiplication of two paths consists of joining them beginning to 
end to form a new path(Fig.2). The U's can then be thought of as representations 
of the abstract group relations given by eqn's (2. 4). It is easy see that, on the 
lattice, (2.4b) is trivially satisfied if we build up the path(string)link by link; on 
the other, the cancellation condition (2.4a) is not automatic and remains to be 
implemented. 

Before turning to the problem of non-bactracking random walk, as a pre- 
liminary exercise,we will briefly consider the simpler and extensively studied 
problem of random walk without constraints(13). This problem will be solved 
first by the operator and then by the path integral method. The reason for 
going through this simple exercise is that the methods used generalize easily to 
the real problem. We start with a regular lattice with spacing d in D dimen- 
sions,and label the lattice sites by the vector x,whose components are given by 
Xa- Boldface letters will always denote vectors in the target space. Let the or- 
thonormal set of states |x > be discrete eigenstates of x with the corresponding 
eigenvalue, and let Pa = —ijr- be the components of the conjugate momentum 
p. The operator exp{idjpa) shifts by one link in the positive "a" direction for 
j = +1 and by one link in the negative "a" direction for j = —1. We now 
consider a sum over all paths starting at x and ending at x', each weighed by 
a Boltzmann factor exp{—al), where l=nd is the length of the path. It is given 
by the following operator expression: 

G(x,x') =<x'|(l-2e-"'^^cos(rfpJ)-^|x> . (2.5) 



We are also interested in the continuum limit d ^ 0. This limit is taken by 
replacing the cosine by 

Y,cos{dpa)^D-y^Y.Pl (2-6a) 

a a 

and scaling G by 

G = d^-^e-'"^G, (2.66) 

while keeping /t^ = (e""^ — 2D)/d^ fixed as rf -^ 0. The Green's function G then 
satisfies the Helmholtz equation: 

(«:2-V.Mx,x')=5^(x-x'), (2.7) 

which is a well known result. Here,v = ~ ^aP^a — ~P^- 

The operator formulation of the random walk problem does not easily gen- 
eralize to strings; instead, it is better to use the path integral approach. It is 
convenient to introduce link vectors e^ of unit length in the positive "a" direc- 
tion and vectors e^j = je^ representing unit steps in + or — "a" direction,where 
j = ±1. The correlation function G is then given by 

oo 

G=Y^ ^M^ (2.8a) 

m=0 

where 

/m— 1 m 

J] rf^X(m') n d'^Y'{m")exp{iI^"'^ - (m - l)ad), (2.86) 

m'=l m"=l 

and, 

m 

expiil^"^^) =Y,exp(i ^ P(m').(X(m') - X(m' - 1) - rfe„j(m')). (2.8c) 

e m'=l 

In this equation, integration over P enforces the constraint that adjacent X's on 
the lattice differ by one of the link vectors rfCa j.It is also understood that the 
initial and final points of the path,X(0) and X(m), are fixed at X and X'. 

In the continuum limit, the discrete parameter m is replaced by the con- 
tinuous parameter r = md. Again expanding to second order in d according 
to 

Y^exp{-idearV) ^ exp{ln{2D) - ^P'), (2.9a) 



we have the continuum action 

il= f dT{2iX'.P - p2 - K,^), (2.96) 

which descibes a simple Gaussian model, in agreement with eq.(2.7). Here, k'^ 
is a constant. 

Next, we consider a non-backtracking random walk on the lattice. This is 
a Markov process, and it can be handled by slightly generalizing the operator 
approach we have used for the unrestricted random walk. Now we pick a point X 
on the path,with the previous step leading to it given by eaj and the following 
step by e^^fe (Fig. 3). If the path is non-backtracking, a = b,k = —j is not 
allowed. To take this into account, the state representing X should have the 
memory of the previous step. This is accomplished by attaching the indices "a" 
and "j" of the previous step to the state |X > by defining 

\X,a,j >=\X> ®\a,j >, (2.10) 

where \a,j > are orthonormal states. The two point correlation function (see 
eq.(2.5)) is then given by 

G(X, a, J, X', 6, k) =<X.',b,k\{l- e-^'^T)-^ |X, a,j >, (2.11a) 

where the 2D by 2D matrix T is given by 

Tbk,aj = exp{idpbk){l - 6a,bSk-j), (2.116) 

and Pa = —i-iT- as before. It is easy to check that T does exactly what it is 

OX a 

supposed to do; it moves X by one step in all directions except the backtracking 
one. 

To gain some insight into what is going on, it is useful to diagonalize T for 
fixed p. Defining oj = 2J2aCos{dpa), the four distinct eigenvalues are given by 

A« = | + i(c.2-8D + 4)V2, 



A(2) = ---(c.^- 
2 2^ 


-8D + 4)1/2 


A(3) = + 1, 




A(^) = - 1. 





(2.12) 



The corresponding eigenvectors can be written in the following form: 

\i>=Y,aajW,j >, (2.13a) 

where the a's are constants, and the index i = 1,2,3,4 labels the eigenvectors 
corresponding to the eigenvalues A'-^'^'^'^-'. The equations satisfied by the a's are 

c^af^ = 1 - X''^''^^exp{idjpa), 

T.cSf =0, 



(3,4) 
Y ■ 
a,] 



aa,-i= - exp{-idpa)oial, 

Oi^a^-i = exp{-idpa)al^}. (2.136) 

Since T is not herniitian,the eigenvectors are not orthogonal. This apparent lack 
of symmetry between the past and the future is due to the way we chose to label 
the states, and it does not have any physical consequence. For example, physical 
quantities like the eigenvalues are real. We also note that the eigenvectors corre- 
sponding to i = 1, 2 are non-degenerate, whereas the eigenvectors corresponding 
to i = 3, 4 are each D — 1 fold degenerate. It is useful to expand A*^^'^) to second 
order in d in anticipation of the continuum limit: 

Ad) ^2D-1- ^^^ci2p2, 
2D -2 ^ ' 

A(2) ^ 1 + -J^d^p\ (2.14) 

2D-2 ^ ^ ' 

In the limit d — > 0,the eigenvector corresponding to the eigenvalue A*^^) dom- 
inates over the rest; and therefore, in the continuum limit, we could drop all 
the eigenvectors except the first one. However, we shall see later that this is no 
longer true in the case of strings; accordingly, we postpone the discussion of the 
continuum limit to the next section. In passing, we note that the zeroth order 
term in the expansion of A*^^-* is reduced from 2D to 2D — 1 due to the fact that 
the step in the backtracking direction is no longer allowed. 

Let us now recast eq.(2.11) into the path integral language, with the view 
of generalizing it to strings later. In the unconstrained random walk problem, 
only X and P were needed as variables of integration. However, since we now 



have discrete indices a and j to deal with, additional variables of integration have 
to be introduced. These variables could be either bosonic or fermionic, and in 
one and two dimensional spaces, they are known to be equivalent. (Here we are 
referring to the dimensions of the world line and the world sheet). We find that 
fermionic variables have several advantages over the bosonic ones; for one thing, 
the formulas look much simpler, and also finitely indexed fermions naturally 
span a finite dimensional space due to the exclusion principle. Let us therefore 
introduce fermionic c- number variables fjajim) and riaj{m) to keep track of the 
indices a and j. Then the J*^™^ of eq.(2.8b) can in the non-backtracking case be 
written as 

m— 1 

i/M = J2 (iPK).(X(m') - X(m' - 1)) 

m'=l 

+ 2 Y.iVa,j{m)r]aj{m' + 1) - f]a,j{m' + l)r]a,j{m')) 

+ E Vb,k{^')Mbk,aj{m')vaA^'))- (2-15) 

aj,bk 

To establish connection with the operator formulation of the same problem and 
to find M, the transfer matrix for the rj's should be identified with the matrix T 
of (2.11b). Remembering that Tfefc^aj(m) maps the state \a,j > at step m into 
the state 1 6, /c > at step m+1, we must have 

ri{m + 1) = T{m)ri{m), (2.16a) 

where we have used the matrix notation. Comparing this with the equations of 
motion for rj that follow from the lagrangian of eq.(2.15), 

-{r]{m + 1) - vim - 1)) + M{m)r]{m) = 0, (2.166) 

gives us the relation between M and T: 

-(T(m) -T-^(m-l)) + M(m) = 0. (2.16c) 

Parenthetically, the expression for T given in (2.11b) should be slightly modified: 
The operator pi, should be replaced by the integration variable Ph{m). This is 
why T acquires an apparent m dependence. 

9 



The action of eq. (2. 15), with M given by (2.16c), is still not completely 
equivalent to the operator approach; it has to be supplemented by an important 
boundary condition. The transfer matrix T(m) maps a single fermion at step 
m, represented by ri{m),into again a single fermion at step m + 1. On the other 
hand, once we have a functional integral, the resulting field theory will have all 
the states containing arbitrary number of fermions, up to the maximum allowed 
by fermi statistics. These extra states are unwanted; for example, in the zero 
fermion sector (vacuum), the T matrix is trivial (unity), and this sector therefore 
describes the unconstrained random walk. To eliminate these unwanted sectors, 
we note that the fermion number is conserved, and sectors with different fermion 
numbers do not mix. We can therefore impose the boundary condition that the 
initial state has fermion number one; this condition will be preserved by the 
dynamics. This initial condition is the extra ingredient that supplements the 
action given by eq.(2.15). 

3. Non-Backtracking Strings 

In this section, we will try to extend the results obtained so far to non- 
backtracking strings. We have the same D dimensional lattice of the last sec- 
tion, and as a preliminary exercise, we first start with the ordinary unconstrained 
sum over surfaces traced by a string(string world sheet). At this point, we have 
to define what we mean by a surface. It is easy to think of several reasonable 
definitions; most of these, possibly even all of them, will have the same contin- 
uum limit. Here we choose a simple definition: Consider the coordinate variable 
Xa{Tn,n), parametrized by two integers m and n. The surface is defined by 
requiring that two adjacent points on it are connected by a link vector, hence 
the following equations must hold for some a,j,a',j': 

X(m,n) — X(m — l,n) = de\j{m,n), 

X(m,n) — X(m, n — 1) = de\,,j,{m,n), (3.1) 

where the e's are the same unit vectors defined in the last section, and super- 
scripts (1) and (2) are attached to distinguish between them. These constraints 
can be implemented by introducing the lagrange multipliers P*^^-* and P*^^\ and 
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the partition function Z as a sum over surfaces is then given by 
Z = fH d^X{m, n)rf^p(^) (m, n)d^P^^^ [m, n) 

mn 

y^ expii ^(P''"'^^(m, n).(X(m, n) — X(m — 1, n) — de^Um, n)) 

aja'j' mn 

+ p(2)(m,n).(X(m,n) -X(m,n- 1) -rfe2\,(m,n))). (3.2) 



Here, for simphcity, we have not introduced a term propotional to the area 
(Nambu term) in the exponential, since, as we shall see, such a term is automat- 
ically generated in any case. To arrive at a correlation function anologous to G 
of the random walk problem(see eq.(2.8)), again one has to fix the initial and 
the final X's and sum over the number of steps it takes to get from one to the 
other. Let us now consider the naive(formal) continuum limit of this equation. 
The action in this limit is gotten by expanding in powers of d, and by dropping 
terms that vanish in the limit d ^ 0. It is of course necessary to scale the field 
variables correctly to end up with the correct continuum kinetic energy terms. 
Defining the coordinates on the world sheet by 



and scaling by 



we arrive at the result 



(Ji = md, 0-2 = nd, (3.3a) 

p(i'2) _^ (2D)^/^p(^'^\ 

X -^ {2/Df'^d X, (3.36) 

^I ^Jd'-T. (2.P«(-).^ - iP^H-)n (3.4) 



1=1,2 



A few comments on this result are in order: 

a)Integrating over P reproduces the action for the standard bosonic string, gauge 
fixed in the conformal gauge. However, the gauge fixing terms(ghosts) are missing. 
This presumably comes about because the naive continuum limit misses subtle 
quantum effects. It is easily fixed by the standard procedure of introducing 
the world sheet metric, casting (3.4) into reparametrization invariant form, and 
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finally fixing the gauge (14). 

b)Tlie continuum result is rotationally invariant on the world sheet. This is of 
course essential for being able to generalize it to a reparametrization invariant 
form; since we have a constant background metric on the world sheet, the only 
remnant of reparametrization invariance is a global rotation invariance, and this 
should remain unbroken. There is even a discrete residue of this symmetry in 
the case of the lattice: Eq.(3.2) is invariant under m ^^ n,l ^^ 2. When we 
introduce fermionic variables later, this symmetry will generalize to a symmetry 
between the right and left chiral sectors of the model. 

c)It was pointed out earlier that the definition of a world sheet surface was non- 
unique. In fact, the concept of surface is not central in our approach. All we 
require is that, upon identifying n with a discrete time variable,the coordinate 
vector X(m, n) should trace all the acceptable Wilson lines(loops) for a fixed 
n and variable m. On account of the discrete symmetry discussed above, this 
should be true also for fixed m and variable n. Therefore, in our approach, the 
world sheet can be thought of as a direct product of two sets of Wilson lines, 
traced by X(m, n) for either fixed m and variable n, or for fixed n and variable 
m respectively. 

As it stands, eq.(3.2) does not generate a satisfactory string theory, since 
backtracking Wilson loops are present in the sum. To eliminate them, we pro- 
ceed as in the random walk problem discussed in section 2. To the first set 
of Wilson lines, we attach the fermionic variables r/^ •(m,n), and to the second 

(2) 

set the variables rj^j (m, n) ; these variables "guide" the path so that it is non- 
backtracking. For a single path, the transfer matrix T for the r]'s is given by eq. 
(2.11b). In the present case, we have two transfer matrices T*^^'^^ , each given 
by equation (2.11b), with p^ replaced by P^; (m,n) and P^; {m,n) respectively. 
The partition function is therefore given by 

where 
I = h + l2 + h, 
and. 
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h = ^(pW(m, n).{X{m, n) - X(m - 1, n)) + P^^\mn).{X{m, n) - X(m, n - 1)) 

mn 

h = -^2 (^ ("^' '^)v {m + l,n) — f]^^' (m + 1, 11)1]''^' (m, n) 

+ f]^"^' (m, n)r]^'^' {m, n + 1) — r^*-^-* (m, n + 1)1]^"^' (m, n) j , 
h= Y.Y1 v'^^\m,n)M^^\m,n)r]'^^\m,n). (3.5) 

ni" 1=1,2 

In this equation, we have again adopted the matrix notation and suppressed the 
indices a and j attached to the r^'s.From the definition of T and the equation of 
motion for 77, the following relation between T and M follows(see eq.(2.16b)): 

-(T^^\m,n) - {T^^\m - l,n))"^) + M^^\m,n) = 0, 



2 



t 



(T(2)(m,n) - (T(2)(m,n-1))-^) +M2(m,n) =0, (3.6a) 



2 
where T is given by 

^£3 ("^5 n) = exp{idkPi^''^\m, n))(l - 6a,bh,~j)- (3.66) 



Eqn's (3.5) and (3.6) describe the dynamics of non-bactracking paths in 
the functional approach. As explained at the end of section 2, they have to be 
supplemented by boundary conditions on the fermion number. These boundary 
conditions must ensure that one and only one fermion is attached to each of the 
fixed m and the fixed n trajectories. In the case of random walk, the existence 
of a conserved fermion number was crucial in preserving the boundary condition 
along the trajectory. Here, in contrast, we have two sets of conserved fermion 
number densities N^-^\m, n) and N^'^\m, n), where the first denotes the number 
of fermions propagating along the link connecting (m,n) to (m+l,n), and the 
second those propagating along the link (m,n) to (m,n+l). The conservation 
laws, which follow from (3.5), are 

ArW(m + l,n) =N^^\m,n), N^^\m,n + 1) = N^^\m,n). (3.7) 

The following boundary conditions are then imposed: For a fixed initial value 
rrii of m, we require that 

Ar«(mi,n) = l, (3.8a) 
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and, for a fixed initial value rii of n, the condition is 

N^'^\m,ni) = l. (3.86) 

These conditions are preserved along each trajectory by virtue of the conser- 
vation laws, and they eliminate unwanted sectors, including the sector iV*^^) = 
A^*^^^ = 0, which corresponds to the standard bosonic string of eq.(3.2). 

Is the string theory defined by eq.(3.5) the QCD string? In general, we 
cannot answer this question for a string that lives on a lattice. However, we 
expect to show in a future publication(lO) that in the limit of continuous time, 
still keeping space latticized, the string theory we have constructed is equivalent 
to a gauge theory on a space lattice. This theory differs from the standard 
hamiltonian lattice QCD in having non-minimal interaction terms, although the 
formal continuum limit is the same. It is therefore reasonable to hope that 
the continuous space-time limit of our model, which we will study in the next 
section, is in fact the large N limit of QCD. 

4. The Continuum Limit 

In this section, we shall study the naive(or formal) continuum limit of 
eq.(3.5). The world sheet coordinates (Ti,2 are defined and bosonic variables 
X and P are scaled as before (eqs.(3.3a,b)), whereas the fermionic variables are 
scaled by 

The continuum limits of /i,2 are straightforward to compute: 

ax 






^ dri^^^ 



l''"%y''-k- f^-^' 



1=1,2 



Finding the limit of I3 involves little a bit more work. It is convenient to first 
diagonalize T (see eqs.(2.12) and (2.13)) before taking the d ^ limit, and then 
compute M through eq.(2.16c). In this calculation, only terms up to first order 
in d are needed, since higher order terms do not contribute in this limit. To first 
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order, the eigenvalues Aj of M are 

Ai = -( 2D + 1), 

A2 =0, 

A3 =0, 

A4 =0, (4.3a) 

and the right and the left eigenvectors \i > and < i\ are given as a sum over 
states \a,j > and < a,j\ (see eq.(2.10) for the definition): 

o n _ 1 

|1> = {2Dr'/'Y.i^ + ^d ——P^j)\a,j >, 

■iri 

<1| ={2Dr'''Y.(l + -^P^j)<a,j\. (4.36) 

«i 

The states |2 >, |3 > and |4 > are not needed, since, the corresponding eigen- 
values being 0, they do not appear in the expansion of M. Notice that as T and 
M are not hermitian, the left eigenvector < 1| is not the conjugate of the right 
eigenvector |1 >. Also, the terms T and T^^ in (3.6a) have different arguments, 
and this causes some confusion about what the argument of Pa should be in 
eq.(4.3b). However, replacing m — 1 by m, or n — 1 by n changes the final result 
only by terms of order (P and therefore it is permissible. The matrix M, again 
to first order in d, is given by 

Mi^2 = '-^{l+^d ^-^^Pi'''h + -^Pi'''b). (4.4) 

bk,aj 2D ^ 2D -2 * 2D -2 '^ ^ ^ ^ 

Substituting eq.(4.4) in (3.5) then gives the continuum limit of J3. This term 
can be simplified substantially by defining a new set of rj's by 

(0 = r9n^-i/2v-„(0. 



,(0 =(2D)-V2^^i;j, 



aj 



j 



^(0 = Y^ f ^(0, 



EUaVlHj, (4.5a) 



aj 
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and a new set of rj's through the same equations. Here, / = 1, 2 as before, and 
the index c takes on 2D — 2 values. The real constants fca satisfy the conditions 



-'C,C' } 



/ J Jc,ajc',a "c 

a 

T.fc,a =0. (4.56) 

a 

Eqn.'s (4.5a,b) are a set of orthogonal transformations to a new set of variables 
s*-'-*, f^'-* and w^^\ which diagonalize M, and which give the following simple result 
for I3: 

(4.6) 
Now that we have the continuum limit, let us discuss some of its features: 
a) The theory is rotation invariant in the D dimensional target space, provided 
that we assign the following transformation laws to the fields: Pa^''^\v^''^^ and 
v^^''^^ are D vectors, and the rest, s^'^ and w^^''^\ are all scalars.This assignment 
follows uniquely from the fact that Xa is a D vector. If so desired, one can then 
introduce the Minkowsky metric in the target space and end up with a Lorentz 
invariant theory. 

b)The fields w^J^'^' decouple from the rest of the theory and become free fields. Of 
course, when a metric is introduced on the world sheet to restore reparametriza- 
tion invariance, they will couple to this metric and will contribute to the value 
of the central charge and to the calculation of higher order terms in 1/A^. Since 
we will not study these more subtle issues in this paper, from now on these fields 
will be dropped. 

c)Since the background metric is flat, depending on its signature,we expect to 
have either Lorentz or rotation symmetry on the world sheet. To see how this 
works, we first consider the kinetic energy terms for the field rj in eq.(4.2). From 
the structure of these terms, it is clear that r]^^^ and i]'^'^^ have to be identified 
with the the two chiral (left and right) components of a world sheet spinor. 
Accordingly, we define the following two component spinors: 



V> 



a,3 




la, =(^3'^S)' (4-7a) 
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with a similar definition for s and Va- I2 tlien becomes tlie standard Dirac kinetic 
energy term: 



l2 = tjd^rj:r^a,l>'^, (4.76) 



where /x = 0, 1, and, 



7° = o-^, 7^ = iay, 



r° = (Ti + 0-2, r^ = (Ti - 0-2- (4-7c) 

In this equation, a^^y are Pauh matrices, not to be confused with 0-12, the world 

sheet coordinates. From the structure of the Dirac equation, it is clear that o"i^2 

are the lightcone coordinates, in contrast to r°'^, the Cartesian coordinates. It 

then follows that the world sheet metric is Minkowskian, which we have already 

anticipated by writing the path integral (3.5) in the Feynman form, with a 

factor "i" in front of the action I in the exponential. We also need the covariant 

two-vecor components P^ of momentum in terms of the lightcone components 
pi, 2. 

pO ^ p(l) + p(2)^ pi ^ p(l) _ p(2)_ (4,7^) 

With these assignments of transformation laws for the fields, the second and 
the third terms on the right hand side of eq.(4.6) can be written in world sheet 
Lorentz invariant form: 

al a,n 



E^^^'M'^^i'^ =\Y.-^i'v^Pa,.- (4.^ 



2 

al a,M 

However, the first term on the right hand side of (4.6), proportional to, 

S^\Y.-^^%^\ (4.9) 

is not world sheet Lorentz invariant as it stands. We will return to this problem 

later. 

d)All of the terms but one in the action have dimensionless coupling constants 

and therefore, they are classically scale invariant. The only exception is again 
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S, the Lorentz non-invariant term discussed above. 

e)The action is parity invariant (left-right symmetric) on the world sheet. 

We will now show that the troublesome term S can be cast into a form that is 
both scale and Lorentz invariant. To be able to do this, it is important to realize 
that the form of the action is not unique; one can multiply any one of its terms 
by one of the number density operators constrained to be unity(eqs.(3.8a,b)). 
Expressed in terms of the fermionic variables on the lattice, these constraints 
can be written as 

1 = N^^\m,n) -^ '^f]lj{m,n)rij{m,n) 
= ^(0^(0 + ^ ^(0^(0 + ^ ^(0^(0. (4.10a) 

a c 

In the continuum limit, this equation becomes 

d{s^'\T)s^'\r) +J2^^dHr)v^'\r)) = 1, (4.106) 

a 

where, as explained earlier, the term involving w's was dropped. The factor of 
d comes from the scaling of the fields in passing to the continuum limit. The 
term S, which is quadratic in field variables, can now be converted into a quartic 
term by multiplying it by a suitable combination of terms that appear on the 
right hand side of (4.10b): 

S =i(iv«s(2)s(2)+iv(2)s«s«) 
a 

-* -S7^S {s-f^S + Y.^al^^Va)■ (4.11) 

^ a 

Finally, putting together eqs.(4.7b,d),(4.8) and (4.11), we arrive at the following 
expression for the continuum action: 

I = J d''r(P^ d,Xa + Ci P^ Pa,, + IS Yd,S + IVa Yd.Va 
+ C2 SJ^S S7^S + C3 S7^S VaJf^Va + C^ Va^Va Vb^^^Vb 
+ C5 S^Va Pa,f^ + C6 Va^S Pa,,i + Cj S^^a Val/iS^ , (4.12a) 

where the summation convention for repeated indices is used, and the partials 
are with respect to the variable r. The constants "c" are given by 

ci = C4 = C7 = 0, C2 = C3 = -iAi, 



_ Ai _ Ai(2D-l) 

"'' ~ "2DV2(2Z)-2)' ""' ~ "2DV2(2D-2)' ^ ^ 

We have written eq.(4.12a) deliberately with arbitrary constants c; written 
in this fashion, it is the most general local action invariant under the symme- 
tries listed after eq.(4.6). In the naive continuum limit, these constants have 
the values given in eq.(4.12b). However, the naive continuum limit is clearly 
not the final answer, since under the fiow of the renormalization group, these 
constants will be driven from their initial values given by (4.12b) to some fixed 
point corresponding to different values for these constants. If this fixed point 
represents pure QCD, it is natural to expect it to have the following properties: 
a) It should be conformally invariant. Otherwise, renormalization would intro- 
duce a new scale in addition to the slope parameter which multiplies the action 
and which we have set equal to one. Since we know that pure QCD has only 
a single scale parameter, conformal invariance has to be an essential feature of 
the world sheet field theory. 

b) The constant Ci and at least one of the constants C5 or Ce should be differ- 
ent from zero at the fixed point. The no n- vanishing of Ci is necessary to have 
asymptotically linear Regge trajectories. With ci not zero, we can integrate over 
P and recover the standard quadratic term in X. On the other hand, if both C5 
and Cg vanished, the fermionic and the bosonic sectors of (4.12a) would decouple 
from each other, and the bosonic sector by itself would simply reproduce the 
standard string model with purely linear trajectories. ( Here the labels bosonic 
and fermionic refer to the world sheet variables). Since we already know that 
the standard string model does not describe QCD, it follows that either csor ce 
should be different from zero. If this is indeed the case, the resulting trajectories 
will be in general curved because of the coupling between P and the bilinears in 
fermions. As explained in the introduction, one expects QCD trajectories to be 
curved, and so this is a step in the right direction. 

c)The action given by (4.12) does not look hermitian.Some of this apparent 
non-hermiticity can be eliminated by redefinitions of fields. On the other hand, 
an imaginary value for C2, as given by eq. (4.12b), is a genuine violation of 
hermiticity. However, this is only the initial value of C2 and not its value at the 
fixed point. Since the model should be unitary, we expect C2 to become real at 
the fixed point. A related question of unitarity is whether the ghosts due to 

19 



the Minkowski metric in target space decouple. We have nothing to say about 
this important problem in this paper. In any case, it is clear that to have any 
chance of a viable theory, the generalized Thirring model (4.12a) must have a 
fixed point with the properties listed above. To investigate the fixed points of 
this model is therefore a task of great importance for the future. 

We close this section with a discussion of the continuum limit of the con- 
straints (3.8), which supplement the action of (4.12). We remind the reader that 
these constraints are necessary to eliminate the spurious sectors present in the 
model. They could also play an important role in making the theory unitary by 
eliminating states with negative norm. In order to write them in covariant form, 
it is convenient to introduce the following vector current: 

r^7sYs + Y.^al^^a, (4.13a) 

a 

whose light cone components turn out to be the fermion densities of eq.(3.7): 

iV« ^ j(-) = i(jO-ji), 
iV(2) ^ j(+) = i(jo + ji). (4.135) 

The continuum analogues of eqns.(3.8a,b) are then imposed on j*^^) at fixed 
initial values of ai = r'^+) = |(r° + r-^) and cr2 = t'-^-' = |(t'^ — t^) respectively. 
For these conditions to be preserved along a trajectory, J should be chirally 
conserved: 

d^r = 0, e^^'d^J, = 0, (4.14a) 

or, in terms of the chiral components, 

9+j(-) = 0, 9_JW = 0. (4.146) 

These are then the continuum version of eqn.(3.7). The remaining problem is 
how to impose the analogues of (3.8a,b) on J. In the quantized theory, J is an 
operator, and it is well known that imposing conditions such as eq.(3.8) directly 
on the operators leads to inconsistencies. Instead, one has to impose these 
conditions on the physically acceptable states, using only the positive and zero 
modes of J. If r^"^^ are compactified on an interval of length 27r, we can define 
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the Fourier modes of J by 

j(±){r) = Y,exp{inT^^^^)j\^l (4.15) 

n 

and write down the conditions that a physical state | > must satisfy: 

'^Sl>=0, ^(yi>=0, (4.16a) 

for n > 1 and, 

jfoVl >= I >, 4^1 >= I >, (4.166) 

for n = 0. 

The conservation laws given by eq.(4.14) are the consequence of the invari- 
ance of the action (4.12a) under both vector and axial U(l) transformations. 
Under vector U(l), the fields transform as 

s -^ exp{i6v)s, s —^ s exp{—i6v), 

Va -^ exp{i6y)va, Va "^ Va exp{-i6v) , (4.17a) 

and under the axial U(l), the transformation is 

s -> exp{i-f36A)s, s ^ s exp^i-fsdA), 

Va -^ exp{i-f36A)vA, Va -> Va exp^i^fsSA), (4.176) 

where 73 = i'-f^'j^. J^ and its dual are then the conserved currents resulting from 
these symmetries. We should therefore add these to the list of symmetries of the 
action listed after eq.(4.6). However, these symmetries are not automatic, since 
it is well known that axial transformations of the type given by eq.(4.17b) are in 
general anomalous. Such an anomaly would destroy the chiral conservation laws 
and would prevent us from imposing the constraints given by (4.16) consistently. 
To see how this anomaly comes about, we define a fermionic supermultiplet Z 
built out D + 1 doublets s and Va'- 

Z=il\ (4.18a) 

and using auxilliary fields Wj,;^, we rewrite eq. (4.12a) as 

I = j (fr[P^ d^Xa + Ci Pii Pa,^ + C5 S^Va Pa,>. + Cg Va^S Pa,^ 

+ iz rd^z + Y^i^Ka^z W,^^ + Wi" IV,,^)), (4.186) 
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where Ki are suitable D + 1 ® D + 1 matrices so that upon integrating over 
the W's, (4.12a) is reproduced. For fixed W, there is an axial U(l) anomaly 
proportional to the trace of the field strength built out of W. The condition for 
the absence of this anomaly is that 

tr{Ki) = (4.19a) 

for all" i".At the fixed point, this implies the following two equations between 
the couphng constants that appear in eq.(4.12a): 

C2 = D^a, C3 = -2Dc4. (4.196) 

Clearly, the argument that led to (4.19b) is only semiclassical and may be mod- 
ified by higher quantum corrections. We went through this exercise merely to 
demonstrate that chiral U(l) invariance imposes further conditions on the pa- 
rameters of the model. 

5. Conclusions 

The main result of this paper is the introduction of a new string model. 
This new model was derived by restricting the world sheet to surfaces traced 
by non-backtracking paths on a lattice and then taking the continuum limit. 
It was argued that such a restriction was the necessary minimum condition for 
the string theory to represent a non-abelian gauge theory in the large N limit. 
To implement this restriction on the allowed paths, fermionic variables on the 
world sheet were introduced, and it was shown that the world sheet action of 
the resulting string theory is the generalized Thirring model of eq.(4.12a). This 
action has to be supplemented by the conditions on admissible states given 
by eqs.(4.16a,b). It was also argued that the coupling constants of the model 
should correspond to a fixed point, and this fixed point should satisfy a number 
of conditions, such as conformal invariance, in order to make connection with 
QCD. 

Much remains to be done to find out whether the model is viable. First, one 
has to establish the connection between the latticized version of the string model 
and QCD on the lattice(lO), and then, to see whether the continuum string 
makes sense, one has to find out whether the generalized Thirring model(4.12a) 
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has a fixed point with the required properties. Hopefully, the work presented 
here will stimulate research along these directions. 
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